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ABSTRACT 

Structural alloys embody internal mechanisms that allow recovery of 
state with varying stress and elevated temperature, i .e. , they can return to 
a softer state following periods of hardening . Such material behavior is 
known to strongly influence Structural response under some important ther- 
momechanical loadings, for example that involving thermal ratchetting. 

Here, we investigate the influence of dynamic and thermal recovery on the 
creep buckling of a column under variable loading. The column is taken as 
the. ideal ized (Shanley) sandwich column. The constitutive model, unlike the 
commonly employed Norton creep model, incorporates a representation of both 
dyn^ic and thermal (state) recovery. The material parameters of the con- 
stitutive model are chosen to characterize Narloy Z, a representative copper 
alloy used in thrust nozzle liners of reusable rocket engines. Variable 
loading histories include rapid cyclic unloading/reloading sequences and in- 
termittent reductions of load for extended periods of time; these are super- 
imposed on a constant load. The calculated results show that state recovery 
significantly affects creep buckling under variable loading. Failure to 
account for state recovery in the constitutive relations can lead to noncon- 
servative predictions of the critical creep buckling time. 
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INTRODUCTION 

The influence of dynamic and thermal recovery on the high temperature 
behavior of structural alloys is wel 1 , recognized [1,2] (see Figs. 1 and 2). 
For example, recovery is believed to play a major role in one of the central 
structural problems relating to liquid nietal breeder reactor design, the 
problem of thermal ratchetting. Failure to account for recovery effects in 
structufaV analyses involving repeated thermal transients has been shown to 
give cjualitatively incorrect and rioncbnservative predictions of ratchetting 
in some instances [3]. Special provisions for taking recovery into account 
(e.g. , the so cal led d- reset procedure) [4] a(re now finding their way into 
documents for guiding structural analysis .in the U.S. nuclear industry. The 
effects of thermal ratchetting have also ’been observed in thrust nozzle 
liners of reusable rocket engines. ’ 

A second area in which recovery effects are thought to play a primary 
role, and one that impacts essentially all high temperature system design, 
is that of creep crack growth under creep/fatigue (variable stress) condi- 
tions. Creep (state) recovery is believed to be one of the main causes of 
acceleration of crack growth in the creep regime and under variable loads. 

A study by Kubo [5], using the same constitutive relations [2,6] employed in 
the present work, has shown that the experimentally observed acceleration of 
creep cracks under variable stress can be predicted through the inclusion of 
appropriate mechanisms of recovery in the constitutive relations. 

In this paper we discuss a third area of mechanics in which recovery is 
expected to have significant influence, that is the problem of creep buck- 
ling under variable loading. The presence of mechanisms of recovery 
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can allow creep strain rates to increase following periods of hardening and 
thus cause acceleration of creep buckling. Load reversals may be accompanied 
by dynamic recovery as in Fig. 1 and load reductions by thermal recovery as 
in Fig.' 2; in either case, creep rates are increased upon reloading, tending 
to reduce the critical time to buckling. 

In the vast majority of creep buckling studies the constitutive model 
used has been of the Norton type [7-9] in which the (steady state) creep 
rate is taken as a function of stress and temperature alone. Some investi- 
gations have made use of hardening theories [10-12] that include a primary 
creep phase. However, these theories do not generally allow for inelastic 
state recovery and, consequently, for rejuvenation of primary creep. These 
constitutive theories are adequate in the case of constant loading but may 
significantly over predict the time to buckling under certain types of vari- 
able loading. 

Here, we examine the elevated temperature creep buckling problem under 
variable loading using the simple Shanley model [13] of a column but with a 
constitutive model [2,6] that embodies a representation of both dynamic and 
thermal recovery. The constitutive model is that developed principally by 
Robinson at Oak Ridge National Laboratory and NASA Lewis Research Center 
and, as indicated earl ier, was that used in the studies concerning creep 
crack growth by Kubo. 

The present study was motivated partly by the occurence of what appears 
to be local ratchetting-buckl ing failures in the throat liner of the main 
thrust nozzles of reusable rocket engines, notably the NASA Space Shuttle 
main engine. For that reason the material parameters used in the constitu- 
tive equations are chosen to represent a copper alloy, Narloy - Z, which is 
typical of materials used in rocket engine thrust nozzle liners. 
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We shall first state the constitutive model and, secondly, specify the 
geometry of the Shanley column model- We then investigate creep buckl ing 
behavior under constant loading; cyclic loading which includes relatively 
rapid load reductions and reversals (dynamic recovery); and cyclic loading 
which includes intermittent reductions of load for extended periods of time 
(thermal recovery). Some of the calculations are repeated fqr comparison 
using a roughly equivalent Norton type creep model. A discussion is given 
of some physicaT aspects associated with the constitutive model and the re- 
sults and finally, we state the conclusions drawn from the study. 


THE CONSTITUTIVE THEORY 

As indicated earlier we make use of the constitutive law reported ip 
[2,6] where a complete multiaxial statement of the model is given. Here, 
we state the model in uniaxial terms together with the corresponding para- 
meter values appropriate for the copper alloy Narloy-Z. 


Flow Taw: 



n 

AF sgn(o-a) ; F > 0 and o(o-a) > 0 
0 ;F<0orF>0 and o(o-a)<0 


( 1 ) 


Evolutionary law: 




^ ■ .n , 

B - R|a| Sgna; 3nd oa>0 


:p 


m-B 


B e - R a Sana; a<a or oa<0 

“o 0 y ’ 0 


( 2 ) 
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(3) 


in which 


F = - 1 


1.60 X 10-8 

4 

8.73 

1.07 X 10-8 (y2) + 1 60 

(4) 

209.6 - 0.20 (T) 

1.46 X 10 7 

1.06 X 10-7 

0.2 

•p 1 

e indicates the inelastic strain rate in h"-*- and o the applied uniax- 
ial stress in MPa. In the general form of the theory, < (MPa) is a Bingham^ 
Prager threshold stress playing the role of a scalar state variable, here it 
is taken to have the constant value given above. For present purposes, 
a (MPa) is the single inelastic state variable . It represents the uniaxial 
component of a tensorial state variable (internal stress) that appears in 
the multiaxial formulation of the theory, og is the minimum attainable 
value of |ci| and plays a primary role in the representation of dynamic re- 
covery (see Di scussion) . The remaining parameters, in some cases functions 
of temperature (T), are consistent with the units MPa, h and degrees Kelvin. 
The elastic response is characterized by a linearly temperature dependent 
Young's modulus E (MPa) given by 


and A = 

n = 
m = 
S = 
k2 = 
H = 
R = 
“o = 


E = 1.47 X 108 - 70.5 (T) 


(5) 



Equations (1-4) incorporate both dynamic and thermal recovery and therefore 
are capable of predicting the type of behavior illustrated in Figs. (1) and 
( 2 ) . 

The parameter values for Narloy - Z were determined from uniaxial ten- 
sile and stress relaxation test data. As the complete data base normally 
required to characterize a particular alloy using the present constitutive 
theory was not available the tensile and relaxation data had to be consider- 
ed sufficient. 

Tensile data were found [14] over a wide temperature range (~30-81lK) 
but at only a single strain rate 0.002/s. Limited stress relaxation data 
[15] were found at the temperature 811K (lOOOF); these were used in the 
absence of creep data. The most serious data deficiency in the present 
context concerns the lack of tests giving a direct measurement of recovery 
effects, i.e., tests such as strain or stress transient dip tests or open 
loop cyclic tests involving partial stress reversals. Nevertheless, it is 
believed that the material characterization is adequate for a reasonably 
quantitative study of creep buckling behavior under variable stress. 

Figure 3 shows a comparison of the tensile data with predictions based 
on equations (1) - (4). Figure 4 gives a similar comparison of typical re- 
laxation responses with a prediction. Neither figure is intended as a 
demonstration of the predictive capability of the constitutive model but 
simply as an ass'essment of the correlation with existing data. The predic- 
tive capability of the constitutive theory has been adequately demonstrated 
relative to other alloys in earlier publications [2,6]. 
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THE SHANLEY COLUMN 


The column model adopted here is the sandwich idealization introduced 
by Shanley [13] and used by Kachanov [16] and others. The column geometry 
is indicated in Fig. 5. All of the deformation is presumed to occur in the 
slender bar elements 1 and 2 of length h and cross sectional area A/2. The 
remainder of the column of length L (L>>h) remains rigid. The width of the 
column, i.e., the distance separating bars 1 and 2, is also taken as h. 

The column is loaded by a time dependent load P(t) (shown positive); 
the lateral displacement of the load point at any time is denoted by u(t). 
The stress and total strain in bars 1 and 2 are denoted by oj, 02 and 
ei = S^/h and 62 = ^ 2 /^’ '"espectively. 

Equilibrium requires that 


o + 0 = 2o« 

2 1 ° 

and 


( 6 ) 


4a u 



where Oq = P/A. 

Differentiation with respect to time provides the rate form of the 
equilibrium equations, 


and 


0+0 = 2 on 

2 1 ° 



4u • 
h °o 


+ 40jj 


h 


( 8 ) 

(9) 


From geometrical considerations (Fig. 5) the compatibility condition 
is obtained as 
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( 10 ) 


( 11 ) 


in which Uq is the initial displacement (imperfection). The rate form of 
the compatibility relation, is 



Deconpbsing the total strain rates in each bar into elastic and inelas 
tic contributions gives 




(13) 



(14) 


where E is the Young's modulus as specified in equation (5). The inelastic 

• p -p 

strain rates and are, of course, obtained by applying the inelastic 
constitutive equations (1) ^ (4) to each bar. 

Combining the equilibrium equations (8) and (9), the compatibility 
equation (12), and the constitutive relationships (13) and (14) leads to 


. 4L . ' M j. 4L 2u, 2L ,'P ;p v 

0 ^ (1 + -FTT 0 ) = O (1 + rT- o - — r-) - 0 ( E, - ) 

1' Eh o' o' Eh 0 h' h o' 1 2 ' 


(15) 
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Under constant compressive stress (i.e., Og <0 and Oq = 0) , and in 

•p *p 

the absence of inelasticity (e^ = =0), equation (15) gives the Euler 

critical stress for the idealized column as 


Further, calling 


Eh 

’c ■ 4L 


(16) 


p = 


(17) 


and 


n - 2u 

- -F 


(18) 


where n is termed the nondimensional displacement, we have from equation (15) 


i ^;p _ ;p 


°1 = 


(19) 


Now, from equation (6) 


®2 " 2og - 


( 20 ) 


and from equations (12), (13), (14) and (18) 


" = ll- ( r - ^2 ) 


( 21 ) 


The coupled system of equations (1), (2), (3), (19), (20) and (21) to- 
gether with the appropriate initial conditions, including the initial imper- 
fection Hg = n(0), allow the (nondimensional) displacement n(t) to be calcu- 
lated for a specified temperature T, Euler critical stress og and history of 
loading Oo(t). The results of several such calculations for various loading 
histories are presented in the following section. 
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RESULTS 


rAll of the calculated results are isothermal with the temperature taken 
to be 811K (lOOOF). The Euler critical stress is 200 MPa and the loading/ 
unloading ramp rate |oq| is 96 MPa/sec. The governing system of equations 
was integrated using a self-adapting Adams-Bashforth predictor-corrector 
method with a 4th order Runge-Kutta method as a starter. The calculations 
were performed in double precision on a Prime 850 computer with an upper er- 
ror bound of 10"^. 

In all cases, the criterion defining the critical time to buckling t^ is 
taken as 



( 22 ) 


This is equivalent, under constant load, to the criterion 


0=0 

1 


i.e., the stress in bar element 1 becoming zero. 


( 23 ) 


Behavior Under Constant Load 

In all calculations the compressive load is first applied from zero, 
where the deforming elements are considered in a virgin state (i.e., 
aj_ = Q2 “ ® nominal value of the applied stress Oq = P/A of -35 MPa. 

Thus,. with the Euler stress oq = 200 MPa the nominal value of the ratio p 
(equation 17) is 



-0.175 


( 24 ) 
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In this section we present the calculated results for a constant load, 
i.e., with p held constant at the above value. Figure 6 shows the response 
n(T) under these conditions for four different values of initial (nondimen- 
sional) displacement no = 0.02, 0.05, 0.10, and 0.2. The time t = t/tg in 
Fig. 6 is nondimensional, being normalized with respect to the critical time 
corresponding to Hq = 0.02, i.e., to “ 6h. 

It is seen in Fig. 6 that the critical time to buckling is reduced by a 
factor of almost 8 with an order of magnitude increase in the initial dis- 
placement. 

For the sake of comparison, all subsequent calculations are taken to 
have the initial imperfection no = 0.02 and are presented in terms of the 
nondimensional time t = t/to- 

Behavior Under Variable Load 

We first consider the effect of rapid load reductions and reversals 
superimposed on the constant load Oq = -35 MPa (p = -0.175). As shown in the 
inserts of Fig. 7, load interuptions occur at time intervals of to = 0.15 
with varying amplitude. These include reductions (in the tensile direction) 
of Oo to -20 MPa and to 0 MPa (insert a), to +20 MPa (insert b), and a com- 
plete reversal to +35 MPa (insert c). The calculated effects of these his- 
tories on the creep buckling response n(t) are shown in the respective curves 
a, b, and c of Fig. 7. 

Curve a corresponds to the loading histories indicated in insert a, i.e. 
with reductions to -20 MPa and 0 MPa. The response curve for each is identical 
to that of Fig. 6 for the same initial displacement 0.02, indicating that no 
change in the buckling response has occurred. Response curve b corresponds 
to insert b in which there is a load reversal to Oq = +20 MPa. Here, we begin 
to see a measurable change in the calculated creep buckling response, the 
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critical time having diminished from the constant load case by about 10%. 
Finally, in curve c, we observe a reduction of more than 30% in the time to 
buckling. This corresponds to the history of insert c where a complete load 
reversal from Oq = -35 MPa to Oq = +35 MPa occurs at each load interruption. 
Immediately following each reversal, we see evidence of the reappearance of 
primary creep. This is attributed to the presence of dynamic recovery (Fig. 
1) where creep (or relaxation) is observed to be accelerated with stress re- 
versals - even in the absence of significant reversed inelastic strain. Al- 
though, to the knowledge of the authors, these effects have not been observed 
directly in creep buckling phenomena, it is expected that such effects can 
occur, on the basis of the experimental observations illustrated schematical- 
ly in Fig. 1. 

The influence of dynamic recovery is best understood by considering the 
state space (o,a) of Fig. 8. Note that in Fig. 8 and all subsequent repre- 
sentations of the state space, compressive o and a are shown as positive and 
are plotted upward and to the right, respectively. Hereafter, the relevant 
quadrants of the state space will be referred to as the first (oa>o) and the 
fourth (oa<0). In Fig. 8 the trajectory of the state point ( 02 , 02 ) 
bar element 2 is shown corresponding to the constant load response curve a 
in Fig. 6, The segment OA traces the path of the state point during initial 
load-up to Oq = -35 MPa. Some inelasticity is indicated over path OA by the 
increase in the inelastic state variable 02 . As Oq is held constant, bar 2 
creeps under nearly constant stress and the state point moves toward B. As 
the geometric nonlinearity becomes prevalent the stress in bar 2 increases 
(as that in bar 1 decreases) and the state point moves toward C. Point C 
corresponds to the buckled condition n = 1 in Fig. 6 (curve a). 
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Similarly, Fig. 9 shows the trajectory of the state point for the load- 
ing; histories illustrated in insert a of Fig. 7, i.e., involving load reduc- 
tions to Oq = -20 MPa and to oq = 0. Here, we see the effect of the abrupt 
load changes as vertical (elastic) trajectories in the state space, resulting 
in no overall change in the state path OABC from that just considered for a 
constant load. This, of course, results in the same creep buckl ing response 
observed earlier for the constant load. 

The state path of Fig. 10 relates to the loading history of insert b in 
Fig,. 7. Here, we begin to see evidence of state recovery . The load varia- 
tions now produce stress reversals in bar 2 and the state point trajectories 
are not just vertical (elastic) lines as before, but now follow curved paths 
as the state point penetrates into the fourth quadrant (oa<o) corresponding 
to a reversal of stress. The state recovers with the stress reversal, re- 
sulting in a relatively softer state (smaller 03 ) upon reloading. Corres- 
p.pn^ingly, the creep rate is increased in response to each load cycle. Evi- 
dence of increased creep rate following reloading is apparent in curve b(Fig. 7) 

Load histories involving larger stress reversals produce increased state 
recovery. Figure 11 shows the state path in bar 2 for the fully reversed 
load depicted in insert c (Fig. 7). Here, the effect of dynamic recovery at 
each unloading/reloading cycle returns the state point well back into the 
primary creep regime. As observed earlier, clear evidence of the rejuvena- 
tjon of primary creep after each loading cycle is seen in curve c, with the 
;tot,al effect of diminishing the critical buckling time by about 30%. 

As the period of application of the rapid loading cycles (tq = 0.15) 
was chosen quite arbitrarily, we now investigate the effect of the frequency 
of load cycles. Figure 12 shows the buckling response curves h(t) for load 
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" histories involving a full reversal, i.e., oq = -35 MPa to oq = +35 MPa (as 
in insert c Fig. 7), with periods to, tq/2, tq/ 4 and tq/ 8 (the constant load 
response curve is also shown for reference). We see the pronounced effect of 
more frequent (shorter period) stress reversals. Reversals with period tq/ 8 
reduce the time to buckling by more than a factor of 3. 

Next, we examine the influence of load reductions of extended duration 
on the creep buckling time. The loading histories considered are shown in 
the inserts of Fig. 13. Case a (i.e., insert a and response curve a) 
amounts to one of those considered previously (insert a of Fig. 7) where 
the time duration at the reduced load Oq = 0 is effectively zero. This his- 
tory, as noted earlier, produces no change in the buckling response over 
the constant load case. 

Insert b depicts the history where the load, having been applied for 
period Tq, is abruptly removed (at a rate Ioq] = 96 MPa/sec) and held at 
zero for a comparable period tq. The sequence is then repeated. This his- 
tory reduces the critical time as shown in curve b, by about 15%. Here, the 

time T includes only that time in which the load Oq = -35 MPa is applied. 

As before, this behavior is best visualized in the state space ( 02 , 

Fig. 14. State recovery is observed as, in time, the state point moves at 
zero stress toward smaller 02 , e.g., on the first load reduction, from point 

D to E. Reloading returns the state point to a softer state than before the 

load reduction and correspondingly to a higher creep rate. Repetition of 
the sequence thus causes acceleration of the creep buckling process. 

Insert c of Fig. 13 shows a loading history in which the time at zero 
stress is now increased a hundredfold to IOOtq. The corresponding response 
curve c shows a further decrease in the critical time, here about 30%. The 
related state point trajectory shown in Fig. 15 is qualitatively similar 
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to that of Fig. 14 but now shows significantly increased recovery with the 
hundredfold increase in hold time at the reduced stress. 

Comparison with Predictions by Norton Law 

In this section we compare the results presented above with those for 
identical loading histories based on a classical creep law of the Norton- 
Bai ley type, i .e. , 

= B sgn (o) (25) 

A roughly equivalent representation was obtained by calculating the steady 
state creep rates for various stress levels using equations (l)-(4) and 
"fitting" equation (25) to these "data" by choosing optimal values of B and 
N in a least squares sense. This process yielded the following 

B = 2.03 X 10-19 
N = 9.375 

*P V 

which are consistent with units of a in MPa and e in h"-^. Using equation 

(25), the critical time to buckling corresponding to the constant load case 

of curve a Fig. 6, turned out to be t' = 7.8 h, slightly greater than the 

0 

earlier reference tg = 6 h. Thus, in the calculated results presented here 

(i.e.. Fig. 16), the nondimensiona! time t is obtained by normalization with 

respect to t'. The time plotted is that for which the load Oq = -35 MPa is 
0 

applied. On this basis, calculations of creep buckling response n(t) for 
all of the loading histories presented earlier have been included in Fig. 

16. As expected, the time to buckling is completely unaffected by any of 
the variable loading histories. This is because the classical Norton-Bai ley 
respresentation (equation (25)) fails to account for state recovery in any 
form. 
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DISCUSSION 


Dynamic recovery, as illustrated in Fig. 1 and by the state point 
trajectories of Figs. 10 and 11, is characterized in the present constitu- 
tive model through the dual analytical forms of the evolutionary law, equa- 
tions (2). This description is consistent with the viewpoint of Onat [17] in 
representing inelastic behavior through the specification of analytically 
different mathematical forms corresponding to various regions of the state 
space / Here, different analytical forms are specified depending on whether 
the state point lies in the first (third) quadrant oa>o or the fourth 
(second) quadrant aa<o of the state space. In effect, this permits the 
state point to rapidly recover upon reversing the stress, even in the ab- 
sence of significant reversed inelastic strain (Fig. 1). This formulation 
is intended as an idealization of a physical process whereby the dislocation 
structure (or the associated structure of internal stress) is abruptly al- 
ter^ed with a reversal of the applied stress, because previously immobilized 
dislocations are remobilized on their slip planes [18]. As the inelastic 
state variable a is taken as an averaged, phenomenological measure of the 
dislocation microstructure (or its associated internal stress state), it too 
should reflect a rapid change as the stress is reversed. 

The specific functional form of the first of equations (2) is made in 
accordance with the experimental results of Mitra and McLean [19]. The se- 
cond equation (2) can be considered an analytical continuation of the first, 
evaluated at Oq (small a) in the first quadrant of the state space, into the 
second quadrant; or with symmetrical response in tension and compression, at 
-Oq in the third quadrant into the fourth quadrant. This amounts to a 
highly idealized representation of the underlying physical process but. 


16 



nevertheless, captures the essential feature that abrupt mi crostructural 
rearrangements occur with stress reversals. 

The path of the state point in the fourth quadrant oo<o in Figs. 11 and 
12, and thus the extent of recovery , depends on the total strain rate and, 
critically, on the parameter oq. Ideally, oq , or more comprehensively a func- 
tion of o and a replacing it, should be determined from experiments in which 
the recovery of creep rate or flow stress under stress reversals is measured 
directly and not inferred indirectly from available monotonic tensile, creep 
and/or relaxation data as was done here. Nevertheless, the present represen- 
tation is believed to be adequate and consistent with the objective of demon- 
strating the strong influence of state recovery in creep buckling. A more 
comprehensive description of dynamic recovery, in the same spirit as that des- 
cribed, and the relevant experimentation are topics of continuing research. 

Thermal recovery, as depicted in Fig. 2 and in the state paths shown in 
Figs, ,14 and 15, is manifest in the second (negative) term of the Bailey- 
Orowan evolutionary equations (2). In the applications considered here the 
state point recovers in time under constant (zero) reduced stress, giving rise 
to an increased creep rate on reapplication of stress (Fig. 2). Physically, 
this macroscopic behavior is associated with thermally activated, diffusion 
controlled microscopic processes such as climb of edge dislocations, allowing 
dislocations, in time, to bypass immobilizing obstacles, thus producing a 
softer state (smaller a). The important material parameters in equations (2) 
are R, m and B. Ideally, these parameters are determined from both creep data 
and information obtained from stress or strain transient dip tests [6] that 
provide a direct measurement of thermal recovery. Again, as these data were 
not readily available for the alloy Narloy-Z, the pertinent parameters were 
inferred indirectly from available data. As discussed earlier, this approach 
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although not optimal, is considered consistent with the present objectives. 

CONCLUSIONS 

We have examined the creep buckling response of an idealized (Shanley) 
column under some special variable loadings. The two types of loading con- 
sidered amount to superpositions of the following load sequences on a con- 
stant applied load: 

1) rapid cyclic unloading/relpading sequences involving stress rever- 
sals, and 

2) cyclic loading that includes intermittent reductions of load for 
extended periods of time (at temperature). 

Although the sandwich column model used is highly idealized, the con- 
stitutive model is quite comprehensive, in as much as it incorporates a re- 
presentation of dynamic and thermal (state) recovery. There is substantial 
experimental evidence that many structural alloys embody internal mechanisms 
at elevated temperature that allow inelastic strain rates to increase (re- 
cover) following periods of hardening. In particular, this is believed to be 
true for the representative copper alloy Narloy-Z characterized here. 

The loading sequences examined are not intended to represent prototypi- 
cal loading histories for any particular structural component; instead, they 
were chosen to best illustrate the generic influence of both dynamic and 
thermal recovery on structural behavior in the presence of a creep induced 
instability. It is expected that qualitatively similar behavior will accom- 
pany more realistic variable loading conditions, (and more realistic struc- 
tures), and that neglect of state recovery effects in such cases will lead 
similarly to nonconservative predictions of the critical time to creep 
buckling. 
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We now list the conclusions drawn from this study: 

• State recovery (dynamic and thermal) can have a significant effect 
on Creep buckling behavior, i.e., on the critical time. 

• Failure to account for state recovery in the constitutive equations 
can lead to nonconservative predictions of the critical buckling 
time under variable loading. 

• A classical Norton-Bailey type creep law as cotimonly used in creep 
buckling analyses, does not account for state recovery and therefore 
may significantly over predict the time to creep buckling under 
variable loading. 

• It is important that constitutive models which are to be used in 
creep buckling analyses involving variable loading allow for recovery 
effects, and furthermore, that the characterization tests used for 
determining the pertinent material parameters include direct measure- 

‘ ments of state recovery. 
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Schematic representation of dynamic recovery. In the 
absence of stress reversalsr the stress relaxations 
over a fixed time ab,cd and ef show successive 
hardening. After the stress reversal fgh, stress 
relaxation hj shows evidence of (creep) softening, i 
evidence of recovery of state (see [1]). 



Schematic representation of thermal recovery. Response 
in an interupted creep test typically exhibits 


relatively small strain recovery 


but measurable 


softening ( depending on tlie interval ) following a 
stress reduction, i.e. recovery of state (see [2]). 


stress (MPa) 


a) b) 



strain (%) strain (%) 

. Fig. 3 Comparison of a) tensile data [14] and, b) predictions 

for Narloy Z at various temperatures. Strain rate is . 002/s. 


24 





(BdN) 


Comparison of stress relaxation data [15] and a prediction 
for Narloy Z at 811K (lOOOF). Cross hatched region Indicates 
range of three tests. Starting stress Is 120 MPa. 
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Fig. 6 Nondimensional displacement versus time ior 
imperfections of H = 0.02 (a), 0.05 (b), 0. 
0.2 (d). ^ 
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Nondiinensional displacement versus time showing tlio 
effect of the frequency of load reversals. Respons 
curves for constant load (-35 MPa) and periods of 1 
to/4 and to/G. 
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Nondimensional displacement versus tine for all loading 
histories using classical Motorn type creep law. 
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